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Abstract 

A recently proposed discrete version of the Schrodinger spectral problem is considered. 
The whole hierarchy of differential-difference nonlinear evolution equations associated to this 
spectral problem is derived. It is shown that a discrete version of the KdV, sine-Gordon and 
Liouville equations are included and that the so called 'inverse' class in the hierarchy is local. 
The whole class of related Darboux and Backhand transformations is also exhibited. 

1 Introduction 

In recent years there has been a growing interest in the field of discrete dynamical systems, i.e. 
systems that can be described by ordinary and\or partial, generally nonlinear, difference equations. 

Such systems arise and play an important role in a very large number of contexts and have an 
extensive range of applications: mathematical physics, chaos, fractals, disordered systems, biology, 
optics, economics, statistical physics, numerical analysis, discrete geometry, cellular automata, 
quantum field theory and so on[]. 

Even if powerful analytic tools were developed in the last decades to deal with difference 
equations, due to the nonlincarity of the systems of interest, no technique of solution is available 
for most of such dynamical systems. Thus it is clear the relevance and the interest of any new 
integrable one. 

In this paper we introduce a whole class of interesting nonlinear differential-difference evolution 
equations which correspond to isospectral deformations of the new Schrodinger discrete spectral 
problem introduced by [Q and investigated in Q and which, therefore, are integrable via the inverse 
scattering method. We show that in this class are included the new discrete integrable version of 
the celebrated KdV found in M, the discrete sine-Gordon and Liouville equations and a whole 
hierarchy of local equations, which as far as we know is new. We give also the explicit recurrence 
operators to construct the Darboux Transformations and the Backlund Transformations for the 
whole class. Of course the Backlund transformations are interesting 'per se' as difference-difference 
dynamical systems and moreover they can be suitable to construct special solitonic solutions for 
the considered differential-difference equations. 

1 The related literature is so large that it is impossible to quote even if in a very partial way. The interested reader 
can find a lot of relevant references in J. of Phys. A, 34 n. 48 (2001), Special issue: Symmetries and integrability 
of difference equations (SIDE IV) 
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2 The spectral problem 



Let us consider the spectral problem 

L(n)ip(n; A) = A?/>(n; A) (2.1) 

with 

L(n) = E 2 + qi^E 1 (2.2) 

where n is a discrete variable (n G Z), A G C is the spectral parameter and _E fc is the 'shift' operator 
defined by 

E k 4>(n) =4>{n + k), k = 0, ±1, ±2, ... (2.3) 

This spectral equation is the discretized version of the Schrodinger equation recently obtained by 
Shabat Q iterating Darboux transformations of the continuous Schrodinger equation and its direct 
and inverse problems were first introduced and studied in 

In the following the 'potential' q and consequently the eigenfunction ip will be considered 
depending as well on the continuous time variable t, while the spectral parameter is considered 
time independent, so that we can apply the inverse scattering method and consider the associated 
difffcrcntial-difference nonlinear equations. 

In the following we shall often use the shorthand notation 

<f>(n) = <j>, 0(n + fc) = fe , fe = ±l,±2,... (2.4) 

We will use also the first order difference operators 



A = E 1 - E° (2.5) 

A = E 1 + E° (2.6) 

and their inverses A -1 , A -1 . Useful formulas for such difference operators and for difference 
equations are reported in Appendix A. 



3 Isospectral hierarchy 

We are looking for nonlinear discrete evolution equations associated with the isospectral deforma- 
tions of the discrete Schrodinger operator L introduced in fl2.2| ). They can be obtained from the 
Lax representation 

L=[L,M] (3.1) 

where M(n,t, E k ) is an opportune shift operator and dot denotes time differentiation. We have 

q(n, ^E 1 = V{n, tjE 1 (3.2) 

where V(n,t) = V(q,qk) is a function depending on q and its shifted values. 

To construct the hierarchy of isospectral nonlinear discrete evolution equations we can use a 
sort of dressing procedure. Precisely, we look for the recursion operators M. and C which allow 
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to construct new admissible operators M' and V in terms of known M and V. Following the 
technique introduced in H and H , we start with the ansatz 



M' = LM + AE° + BE 1 



(3.3) 



where E° is the identity operator and A and B are functions, depending on the q and q k , to be 
determined in such a way that 

V'E 1 = [L, M'] (3.4) 

or explicitly 

V'E 1 =q(A 1 - A) E 1 + (A 2 -A + qB 1 - q x B + qV\ ) E 2 + (B 2 - B + V 2 ) E 3 . (3.5) 

Now imposing for compatibility the vanishing of the terms in E 2 and E 3 and setting, for 
convenience, 



we get the two conditions 



C = A X -A 



B 2 -B = -V 2 

d + C = qi B- qB 1 - qVi 



The general solution of the second order difference equation (|3.7| ) reads 

oo 

B = b + b(-l) n +Y,V 2k+2 



(3.6) 



(3.7) 
(3.8) 



(3.9) 



fe=0 



where b, b are two arbitrary constants. 



Taking into account (3.9) the general solution of the first order difference equation (3.8) reads 



C = C-qJ2 V 2k+1 - 2 J2 (-!)" Ik E Vk + 2 ^ 



k=0 



fe=l 



3=0 



(7 = c(-l) n -6( 9 + 2^(-lf ft j + 6 (-if U + 2 X> 



(3.10) 
(3.11) 



fc=i 



where c is again an arbitrary constant. 



Now, from @, (|3J), ( |3T0| ) and ( |3TT| ) we get 

V = V + cv 

where the recursion operator C is given by 

1 oo oo oo 

<CV = -q\qJ2 V2k+i + 2 £ (-1)" % £ 14+2 J+ i 

J=0 



fe=0 



fc=l 



while 



L + 2^(-lf ft +6(-l) n « + 2^ ft 

V fc=i / V k=i , 



(3.12) 



(3.13) 



(3.14) 
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Note that starting from the trivial operators M = 0, V = we get from (3.12) the non trivial 
starting point V. Indeed, taking into account the fact that b, b, c are arbitrary constants, we have 
three different independent starting points 



V=(-l) n q 

_ / °° N 

V = g U + 2^(-l) fe (Zfc 

V k=i j 

/ oo N 

V=(-l) n q g + 2£> 



(3.15) 
(3.16) 

(3.17) 



so that the class of isospectral nonlinear discrete evolution equations associated with the spectral 



operator (2.1), (2.2) is a superposition of three hierarchies and reads 

q = a{C)V + (3 (C)V + 7 (C)V 



(3.18) 



where a, /3,7 are entire functions of their argument. The possibility of using negative powers of 
the recursion operator will be considered below. 

It is easily seen that the recursion operator C in ( 3.1 3| ) can be written, by using the difference 
operators (2.5), (p. 6]) and their inverses (see Appendix A), in the more compact and elegant form 



C = -qAA~ 1 qE 1 A~ 1 A~ 1 . 
The explicit form of the recurrence operator M. for M (M' — MM 



V 



easily obtained from the above formulas (see (|Q|), (|3.6|), (|3.9p, (3.1C), (3.11)) 



(3.19) 
CV) can be 



Remark 1 Note that all the M operators generated by the recursion \3. 3j ) contain only positive 
powers of the shift operator E. 



Remark 2 Due to the non-local character of the recursion operator C (see (3. IS)) and taking 
into account the starting points (3.15), ( 3.1 dp , (3.1V), all the equations in the class (3. It) are 
apparently non local (except the rather trivial q = V = (— 1) q). However it is easily seen that 



(3.20) 
(3.21) 



q = V = q\q + 2 {~ 1) Qk I and q = V = (—1) q I q + 2 ^ Ik ) imply respectively 
\ k=i J V k=i J 

mi +qq\ = qqi {q-qi) 

<iqi + q<ii = (-i) n 99i (qi + q) 

which are local although not explicit. The first equation was derived and studied in ref. ||/. 



4 The 'inverse' class 



As noted in the previous section, the equations in the class (3.18) correspond to positive shifts 
in the Lax operator M. In this section we investigate the class of isospectral nonlinear discrete 
evolution equations which corresponds to operators M constructed in terms of negative powers of 
the shift operator E, 

These equations can be obtained using the inverse of the recursion operator £, namely 

C' 1 = -E~ 1 AA-A- 1 A-. (4.1) 
q q 
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Thus, given a valid couple M and V, one can obtain a new valid one M' = Ai^M, V' = £~ 1 V. 
Starting from V = and taking into account that an arbitrary constant c is in the her of A we 
have 



V = £ _1 



-2 c 



(4.2) 



1 1 

and the whole 'inverse' class of isospectral nonlinear discrete evolution equations reads 

q = V (£- 1 )V (4.3) 

where rj is an arbitrary entire function of its argument. 

However, in order to study the properties of this class of equations, it is more convenient to 
construct it explicitly through the ansatz 



JV 



m (n) = y2 A (k) E k-N-i_ 



(4.4) 



fc=l 



Inserting the above expression in the Lax equation (3.1) it is easily seen that the functions 
(k = 1, 2, ...N) are determined recursively by 



qA<f +1) -Qs-nA^+V =A^ -4 s) , ^4 (0) =0, s = 0,1,...,N-1. 
The corresponding nonlinear discrete evolution equation reads 



(4.5) 



(4.6) 



The general solution of (4.5) is, according to (A. 14) and (A. 17) 



^(s+l) = 



N-s 



q-j+k+i 



3= 1 



n i-i 

3= 1 



(4.7) 



where the constants a^ s+1 - ) arise from the solution of the homogeneous part of the equation and 
can be chosen equal to zero except — 1 that gives a convenient starting point for the iteration. 

In spite of the infinite series in <K% the are local for any TV. In the Appendix B we 

solve explicitly the recursion relation proving by induction that 



(-1)- 



fci+i 



E E 



E 



fcs-l+l 

E Qfci 



JV-s 

II Q-m kl= - N k 2 =-N+l k j =-N+j-l k 3 =-N+s-l 
m— 1 



(4. 



where 



i 



qkQk+l 

For the reader's convenience, we write here explicitly the first three equations in the hierarchy. 



(4.9) 
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1. For N = 1 we get 

• f 1 1 

q= 

\9i 9-1. 

which was implicitly considered in ||. Note that setting 

1 



9 = 



we get 



2. For N — 2 we get 



U = — U (Ui — M-i) 



1/1 1 \ 1/1 1 



9? V92 5/ 9-i \9-2 9 



which, using the position ( |4.11j) , becomes 

u = u 2 (u\ (li2 + it) — u_ 1 (u_2 + u )) 

3. For AT = 3 we get 



q = (a 2 - a) 

1/1 1 



a 



1 / 1 



OT-i 9-19-2/ 9-i9-2 V99-1 9-19-2 9-29-3 



With the position (4.11) we get the polynomial equation 



u = -m 2 (/3 2 - 13) 

[3 = UU 2 _ 1 (UlU + UU-l + U_l?i_2) + M 2 _ 1 M-2 (lH-1 + U-%U-2 + U-2^-3) 



(4.10) 

(4.11) 
(4.12) 

(4.13) 
(4.14) 

(4.15) 
(4.16) 



(4.17) 
(4.18) 



5 Discrete KdV, sine-Gordon and Liouville equations 



Notice that the first equation (4.10) in the 'inverse' class of the equations related to the discrete 
Schrodinger operator ([2.2]) can be rewritten as 



. , . 9 9i 

9i9 + 9i9 = 

92 9-1 



(5.1) 



and, then, by taking a linear combination with coefficients c and d of this equation with the 
equation ( ft.2C| ) of the direct class we get 



9i9 + 9i 9 = c [ — - — ) + d{qi - q)q 1 q 
\92 9-1.' 



(5.2) 



which is the equation studied in g] and that for c = 2d reduces to a discrete version of the KdV 
equation. 

It was already shown in [pi that the discrete sine-Gordon and Liouville equations are included 
in the hierarchy of integrable equations related to the discrete Schrodinger spectral operator (2.2). 
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They can be recovered by using the special dressing method described above. In fact if, instead of 
starting in the iteration procedure from V — 0, M = 0, we start from 



we get 



M = d t , V 



q = V = -Cq, 



(5.3) 



(5.4) 



which is a tricky equation since in order to extract from it an evolution equation it must be solved 
with respect to q. In order to do this it is convenient to put 



q = -jAAe- 2lp 

with 7 an arbitrary constant and tp — tp(n,t) a new function and 

q = -e 2lfil AP 



with P — P(n,t) to be determined. Inserting (5.5) and (5.6) into (5.4) we have 

pi _ p 2 _ e -2( Vl + V2 ) _ e -2(v>+<pi)_ 

With an opportune choice of the constant of integration and of the sign of P, we get 

P = e^ ( - v+v>1 \ 

Therefore, 

q = -e^-'f - e^-* 2 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



and the differential-difference equation in if can be obtained by imposing the compatibility between 
O) and dpi). We get 



((pi — ip)e Vl v — (tp2 — i Pi)e ipi V2 — je 2ip2 — 7 e 2ip ■ 

Inserting in it 

ipi-<p= - ie -{vi+v) + q> 
with $ = <f>(n) to be determined we derive 

(j) ie -Ol+¥'2) _ $g-(lP+iPl) 

that can be integrated furnishing the evolution equation 



(5.10) 



(5.11) 



(5.12) 



(5.13) 



which, up to a trivial change of function, for 7' = 7 is the discrete sine-Gordon and for 7' = the 
Liouville equation. 

Also the auxiliary spectral problem introduced in |H) for fixing the time evolution of ip can be 
recovered by using our dressing method. In fact we have 



ijj = -M'tp 



(5.14) 
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which, since 

M' = LM + AE° + BE 1 , (5.15) 

can be rewritten as 

i) = -[L, M]i[) + XM^j - Atp - Btp v (5.16) 

Recalling that M = dt and denoting A = — 1 — fc 2 , we have 

k 2 tp = Aip + (B - q)ipi. (5.17) 

Using the formulas obtained for A and B in section |^ and choosing equal to zero the constants of 
integration we get 

A = —A~ 1 qE 1 A~ 1 A~ 1 V (5.18) 

and 

B = -E 2 A- 1 A- 1 V. (5.19) 



Inserting the formulas for q and q obtained in ([k9|) and ( |5.5| ) we recover the auxiliary spectral 
problem derived in S 



k 2 ijj = je-^-^i) - je'^ip!. (5.20) 

The second auxiliary problem introduced in Q is obtained by closing the compatibility conditions 
with the discrete Schrodinger spectral problem (2.1). 



6 Darboux and Backlund transformations 



Let us consider besides (2A) a second spectral problem with a different 'potential', namely let 

= (6.1) 

with 

L(n,t) = E 2 + q(n,t)E 1 . (6.2) 

Let us introduce a Darboux transformation relating ip and ip 

ii^Vijj (6.3) 

where D is an opportune shift operator depending on q and q. This implies a relation between 
q and g, called Backlund transformation, that can be expressed in the following operatorial form 
(see e.g. §, |) 

LV-VL = WE 1 = (6.4) 

where the scalar operator W depends on q and q and their shifted values up to some order. 

Now, following a technique introduced in , we look for the recursion operators Y and which 
allow to construct a valid couple of V and W from the supposed known V, W, that is such that 
V = TV, W = nw. 
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Consider the following ansatz 

V = LV + FE° + GE 1 
where F and G are functions to be determined requiring that 

W'E 1 = LV'-V'L 

for some W. We have 

W'E 1 = {qFi - qF) E 1 + (F 2 - F + qd - q x G + qWi) E 2 + (G 2 - G + W 2 ) E 3 

and imposing the vanishing of the terms in E 2 and E 3 we get the following conditions 

G 2 -G = -W 2 

F 2 - F = qiG -qG-t -qWt. 



The general solution of the second order difference equation (|6.8| ) reads 

G = g + g(-l) n - E 2 A- 1 A~ 1 W 

where g, g are two arbitrary constants. 

The general solution of the second order difference equation (|6.9|) reads 



F = f + f (-1)" + A^A- 1 (q x G - qG x - qWi) 



where /, / are two arbit rary constants. 
Therefore, inserting (3.10) we have 



F = F + F 



where 



and 



F = f + f (-1)" - gE Y A- l A-\q^ -q)- g(-l) n E A~ l £sT x (g_ a + q) 



(6.5) 
(6.6) 
(6.7) 



(6.8) 
(6.9) 



(6.10) 

(6.11) 

(6.12) 
(6.13) 
(6.14) 



F = E 1 A~ 1 A- 1 {-?_! + {q^E 1 - q)E 1 A~ 1 A- 1 }w. 
From @, (|I|), (|t|) and (|6~ll) we get 

W' = W + AW (6.15) 

where 

W=(qFi-qF\ (6.16) 
A = (qE 1 - q) E 1 A~ 1 A- 1 + {q-iE 1 - q)E 1 A- 1 A- 1 Y (6.17) 

Starting from the trivial operators T> = 0, W = we get from ( |6 .15 ) a non trivial Backlund 
transformation W . Indeed, taking into account the arbitrariness of the constants in (6.13) we have 
four different independent starting points so that the class of Backlund transformations reads 



]T a k (A) = 



(6.18) 



k=l 



9 



where the ctk (A) are arbitrary entire functions of their argument and the 'elementary' Backhand 
transformations are given by 



W (1) = {q- q) 

W (2) = (_ 1 )«(gr +(? ) 

W {3) = (qE 1 - q) E 1 A~ 1 A~ 1 - q) 

W {4) = (_!)« (£ E i + q j E^^A- 1 ^ + q). 



(6.19) 
(6.20) 
(6.21) 
(6.22) 



The explicit form of the recurrence operator T for T> can be easily obtained from the above 
formulas. 

Remark 3 In the limit q — > q, as expected, the operator A becomes the operator C 
In fact let us note that A can be rewritten as 



A = (q qE') { £ E*\ qi - qE 1 ) £ E* £ E^qE 1 

3=1 



fc=0 



fe=0 



For q = q we get, recalling (A. It) 

oo 

A 5=9 ={q- qE 1 ) ^(-lfsWA^A- 1 



k=0 



and, thanks to (A.t) we have 



Aq=q = C. 



7 Concluding remarks 

To end this paper we want to outline a number of possible extensions and generalizations. 

First of all we conjecture that the whole hierarchy here introduced is endowed with a double 
Hamiltonian structure. Finding such structure should allow us to exhibit an infinite number of 
commuting conservation laws for the whole hierarchy. 

Moreover, the technique we used to derive our results (recurrence operators, Lax pairs and 
Darboux transformations) can be easily extended to recover differential-difference nonlinear evolu- 
tion equations related with non isospectral deformations of the spectral problem, getting typically 
equations with n-dependent coefficient (see 0). 

The Backlund transformations and the discrete-discrete evolution equations that can be asso- 
ciated to this new spectral problem (see |^|) deserve further investigation as integrable nonlinear 
iterated maps and moreover offer a good starting point for the introduction of new 'integrable' 
cellular automata (see e.g. 0). 

Finally, all these results could be generalized to the non-abelian case considering a matrix 
discrete Schrodinger operator and matrix differential-difference evolution equations (see e.g. [||). 



A Appendix 

For the convenience of readers not so familiar with difference equations and operators we give here 
some useful formulas. 
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Let us define the first order difference operators 

A = E 1 - E° (A.l) 
A = E 1 + E° (A.2) 

where E° is the identity operator. The general solutions of the difference equations 

AX = F, X = X(n), F = F(n) (A.3) 

AY = G Y = Y(n), G = G(n) (A.4) 

can be written as 

X = c + A^F (A.5) 

Y = c(-l)" + A- 1 G (A.6) 

where c is an arbitrary constant, i.e. not depending on n, and the inverses of A and A are chosen 
as follows 



A~ 1 = -J2 Ek (A.7) 

fe=0 
oo 

A- 1 =^2(-l) k E k . (A.8) 

fc=0 

We can also construct higher order operators using these fundamental ones as the second order 
difference operator 

AA = E 2 - E° (A.9) 

and its inverse 

oo 

A^A- 1 = -^E 2k . (A.10) 

fe=0 

Thus the general solution of the second order difference equation 

AAX = F (A.ll) 

can be written as follows 

X = c + c(-l)" + A- 1 A- 1 F (A.12) 

where c and c are arbitrary constants. 

Finally, let us consider the general homogeneous first order difference equation 

AX 1 -BX = 0, X = X(n), A = A(n), B = B(n). (A.13) 

Its general solution can be written as follows 

+00 

F>k 



k=0 



11 



where c is again an arbitrary constant. 

Then for the non homogenous first order difference equation 

AX 1 -BX = F, F = F(n) (A. 14) 

we search a solution of the form 

+00 . 

X = II if (A.15) 

fc=0 k 

with G = G(n) to be determined. We get 

1 + °° R 

a-o-snj' (A16) 

fc=0 



and therefore 



+00 . +00 3 . 

fe=0 K j=0 J fc=0 

B Appendix 

We choose cr 1 ' = 1 and all other a^ s+1 ' = 0. It is convenient to rewrite the recursion relation for 
the y4.( s+1 ) in terms of 

N-s+l 

B W _ A« gf^. (B.l) 



We have 



r( s ) r(«) 

B (-+i) _ = _^L_ + ^ : fi(o) =0 (B.2) 

<79i 9-Af+s-i9-AT+ s 



and then, if we choose all integration constants zero, 



S (s+1) =^o + V ^ . (B.3) 

V 9fc<Zfe+l Qk-N+s-iqk-N+s J 

Now we prove by induction on s that for s > 

-1 fci+i kj-1+1 fc 8 -i+i 

fl (M-D = ( _ 1) . ^ 2 ... 2 ... 2 Q fcl ...Q fcs (B.4) 

fei=-AT/b 2 =-JV+l Hy=— JV+.J-1 fc s =-AT+s-l 

where 

Qk = (B.5) 
qkqk+i 



From (B.3) we have 

= (-1) S+1 A + (-l) s /?2 (B.6) 
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where 



#L = E ^ fel E E "' E "' E Qk 2 +k!+2 ■ ■ ■ Qfe s +fci+2 (B.7) 
fc 1= k 2 = -N k 3 = -N+l k j = -N+j-2 k B = -N+s-2 

oo -1 fc 2 + l fcj-i+l k s -i+l 

02 = ^2 Qki-N+a-1 E ^ ' ' • ^ ••• E Qk 2 +ki • ■ -Qfc a +fei-(B.8) 

fei=0 k 2 = -N k 3 = -N+l k,=-N+j-2 k s = -N+s-2 



In (B.7) we introduce fc' = fci + fej ■ + 2 for j = 2, . . . , s and renaming fc' — » fcj we obtain 



j "'l T n-j r J — ■ ■ ■ ] o aim luuouuiig n,^ — >• n,j 

oo fei + 1 fc s _i + l 

A=E E E E QkrQk 2 ...Qk 3 - (B.9) 

fc 1= fc 2 = -A f +fei+2 k J =k 1 -N+j k s =k ± -N+s 



In (BJ3) we introduce fcj = hi — N + s — 1 and fc^ = fci + fc 3 for j = 2, . . . , s and then rename 



kj — > fcj 

OO fei-S + W fe 2 +l fcj-l + l fcs-1+1 

a= E E E ••• E ••• E 

fci — — iV+s— 1 fc 2 — k\ — s+1 /C3— fci — s+2 fcj — fei +j — s— 1 k s =k\ — 1 

Exchanging the first two sums we obtain 

oo fc 2 +s-l k 2 +l kj-i+1 fc s -i + l 

a= E E E •■• E •■• E 

k 2 = -N fci=max{fc 2 +s-JV,-JV+s-l} fe 3 =fci-s+2 fe 3 =fei+j-s-l fc s =fei-l 

and with fci <-» &2 

oo fei+s-1 fei+1 + 1 fcs-l+1 

E E E ■■• E •■■ E ^ 

k 1 = -N fc 2 =max{fci+s-iV,-iV+s-l} fc 3 =fc 2 -s+2 fe 3 =fe 2 +J-s-l fc s =fc 2 -l 
-1 fei+s-1 fel+1 fcj_l+l fc s _!+l 

= E E E ••• E ••• E ( h ---Q, • 

fei=— iV fc 2 =— W+s-1 fe 3 =fc 2 — s+2 kj=k 2 +j-s-l k s =k 2 -l 
oo fei+s-1 fei+1 + 1 fes-l+1 

+ E E E •■■ E •■■ E Ofa-o*.- (B.10) 

fei— fe 2 — fci+s— AT ^3— fc 2 — s+2 kj=k 2 -\-j— s— 1 fc s =fc 2 — 1 

Let us consider first of all the second term and exchange the sums from the left to the right. After 
j — 1 inversions we have 

OO fe 3 _2+l kj-l+S— (j — 1) fcj_i+l fcj + l+l fc s _i + l 

E - E E E E ••• E <^ 

fci=0 Jy_i=Jfei— JV+(j— 1) fcj=fci+s-A r fcj + i=fcj-s+j fc i + 2 =fcj + (j+2)-s-l fc a =fcj-l 

and exchanging and one S e ts 

oo fej-2+1 kj-i+1 k j + 1 +s-j fcj + i + l fc s _i + l 

E- E EE E ■•■ E ^ 

fei=0 fe 3 _i=fei-iV+(j-l) fe j+ i=fci-AT+j fe^fci+s-jV fc j+2 = fcj + (j+2)-s-l fe s =fe 3 -l 
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so that with kj <-> fcj+i 

oo fej-2+1 fej-1 + 1 fcj+s-J kj+1 k 3 -i+l 

£•■• E EE E ••• E ^ -ft 

fc 1= fe 3 _i=fe!-iV+t3'-l) fej=fci-iV+j fc J + i=fci+s-AT fe i+2 = fe j + 1 + (j + 2)- s -l fe s =fc 3+1 -l 



and we only have to observe that for s = j + 1 the process ends giving for the last sum X)fc S =fc^+ s -jv 



From (B.6), (BJ) and (B.10) it then follows that 



-1 fei+s-l fci + 1 kj-i+1 k 3 -i+l 

B^ = (-iy EE E ••• E ••• E q^-Qk- 

k x =-N fc 2 =— N+s-l k 3 =k 2 -s+2 k j= k 2 +j-s-l k s =k 2 -l 

Again we exchange J2k? an< ^ Sfc 3 getting 

-1 fci+1 k 3 +s-2 fc 3 + l fcj-1 + 1 fe*-i+l 

= (_!)• ^ ^ 53 E ••■ E ••■ E ^ 

k 1 = -N k 3 =-N+l k 2 = -N+s-l fc 4 =fc 2 -s+3 fc i= fc 2 +j-s-l fc s =fe 2 -l 



and for &2 *-* &3 

-1 fei+l fc 2 +s-2 fc 2 + l fcj-i+1 fe s -i+l 

b^x) = (_d- 2 E E ••• E ••• E 

fc 1= — JVfc 2 =— iV+1 fc 3 =-AT+s— 1 & 4 =fe 3 — s+3 fc i= fc 3 +j-s-l fc s =fc 3 -l 



so that after performing all the inversions we recover (B.4) 
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